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Beschreibung 

Verfahren und Vorrichtung zur Lokalisierung von Licht emit- 
tierenden Bereichen 

5 

Die Erfindung betrifft ein Verfahren sowie eine Vorrichtung 
zur Durchfiihrung des Verfahrens zur Lokalisierung von Berei- 
chen, beispielsweise fokalen Lasionen, in einem biologischen 
Gewebeabschnitt, die zumindest wahrend der Untersuchung eine 
10 vom Gewebeabschnitt verschiedene Fluoreszenzeigenschaf t auf- 
weisen, aufgrund derer bei einer Bestrahlung mit Licht einer 
ersten Wellenlange Licht einer anderen Wellenlange emittiert 
wird. 

In Ulnar Mahmood et al., "Near Infrared Optical Imaging of 
Protease Activity for Tumor Detection", Radiology 213:3, 866- 
870 (1999) ist beschrieben, dass f luoreszierende metabolische 
Marker 

• sich entweder ausschliefilich in bestimmten Regionen, bei- 
spielsweise Tumoren, Entzundungen oder anderen bestimmten 
Krankheitsherden, anreichern, oder 

• zwar tiberall im K6rper verteilt sind, aber nur speziell in 
bestimmten Regionen, beispielsweise durch tumorspezif ische 
Enzymaktivitaten und durch zusatzliche Bestrahlung mit- 
tels Licht, aktiviert werden. 

Dieses Prinzip der optischen Fluoreszenz-Bildgebung wird an- 
hand der Figur 1 verdeutlicht , bei dem ein Tumor bei der Be- 
leuchtung mit NIR-Licht (Licht im nahen Inf rarotbereich) 
30 sichtbar wird, nachdem der Maus ein Marker verabreicht wurde, 
dessen Fluoreszenzeigenschaf ten durch bestimmte Enzyme frei- 
geschaltet werden. 

Die Erkennung eines Tumors oder einer anderen markierten Re- 
35 gion geschieht dann liber Bestrahlung der Region mit Licht in 
der speziellen Anregungswellenlange des Fluoreszenzf arbstof- 
fes und Detektion des emittierten Lichtes in der entsprechen- 



15 



20 
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den Emissionswellenlange des Fluorophors. Diese Markerkonnen 
bei Zulassung fiir den Humanbereich beispielsweise bei der 
Krebsf rliherkennung eingesetzt werden. 

Durch Anregen des Farbstoffes mit mindestens einem zeitlich 
veranderlichen Anregungslichtsignal, beispielsweise durch 
zeitliche Variation von Einstrahlungsort und/oder Lichtwel- 
lenlange und/oder Intensitatsmodulation des Anregungslichtes, 
konnen Daten z.B. mittels CCD oder Photomultiplier , Photonen- 
fluss auf der Oberflache des zu untersuchenden Gewebeab- 
schnitts, wie z.B. der weiblichen Brust, an verschiedenen 
Messpunkten zu verschiedenen Variationen gewonnen werden. Auf 
diese Weise erhalt man variations- und ortsabhangige - d.h. 
raumlich zweidimensionale - Messdaten. Im Falle von M Messda- 
ten zu N Variationen sind dies MxN Daten. 

Aus diesen Daten konnen mittels des hier dargestellten Re- 
chenverf ahrens raumlich begrenzte Lasionen - wie z.B. fokale 
f luoreszierend markierte Tumore - mit zum umliegenden Gewebe 
unterschiedlichen Fluoreszenzeigenschaf ten - dreidimensional 
geortet werden. 

Die Methode ist vorgesehen bei Krebs- (Screening-) Untersu- 
chungen der 

• Brust 

• Lymphknoten 

• Schilddruse 

• Prostata 

• intraoperativen Anwendungen 

sowie alien oberf lachennahen Organen, die im Bereich der Ein- 
dringtiefe von Licht liegen und Karzinome (oder andere Krank- 
heiten) entwickeln, fiir die (jetzt oder zu einem zukiinftigen 
Zeitpunkt ) entsprechende Fluor eszenzmarker existieren. 

Zur Fluoreszenzrekonstruktion bzw. Lokalisation existieren 
verschiedene Ansatze. 
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Britton Chance schlug ein Verfahren zur Lokalisation von flu- 
oreszierenden Absorbern in homogenem Medium vor, sogenannte 
Phased Arrays. Dieses Verfahren lokalisiert f luoreszierende 
Inhomogenitaten (Spots) lediglich in absolut homogenen Me- 
5 dien, d.h. Medien mit homogenen Lichtabsorptions- und Streu- 
eigenschaf ten, wie sie in der Anwendung kaum vorkommen, und 
bietet keinerlei Information tiber die Tiefe, in der sich der 
Spot befindet. 

10 Ansonsten wurden verschiedenen Verfahren zur Fluoreszenz- 

Rekonstruktion vorgeschlagen . Bei der Rekonstruktion wird die 
komplette Fluoreszenzaktivitat im gesamten (meist diskreti- 
sierten) Medium ermittelt (ahnlich wie bei nuklearmedizini- 
schen Methoden) , wahrend bei der Lokalisation ausschliefilich 

15 die sich aus dem Hintergrund hervorhebenden Bereiche gesucht 
werden. Rekonstruktionsver f ahren beruhen somit auf der (oft 
iterativen) Losung grofler Gleichungssysteme und sind somit, 
im Gegensatz zur hier vorgeschlagenen, in Echtzeit arbeiten- 
den Lokalisation, sehr zeitaufwendig. Die Rekonstruktionsver- 

20 fahren gehen des weiteren iiberwiegend davon aus, dass das zu 
untersuchende Medium (ahnlich wie bei der Computertomogra- 
phie) von einem Ring von Lichtquellen und Detektoren umgeben 
ist . 

f|5 Einige der bekannten Verfahren sind hier zusammengef asst be- 
schrieben: 

Der US 6,304,771 bzw. der US 5,865,754 sind die Tomographie 
mit f requenzmoduliertem Licht zu entnehmen. Das Rechenverf ah- 
ren benotigt eine Rekonstruktionszeit von 5 Min. auf einem 
30 1GHz Pentium-Rechner bzw. 45 Min. auf einer SUN Sparc 2 Work- 
station. In der WO 02/ 417 60 A2 ist die Tomographie mit Licht 
beschrieben, die ebenf alls eine Rekonstruktionszeit von 5 
Min . auf einem 1GHz Pentium-Rechner benotigt . 

35 Alle diese beschriebenen Verfahren zeichnen sich durch einen 
hohen Rechenauf wand bei relativ kleinen Rekonstruktionsvolu- 
mina aus, eine Berechnung in Echtzeit ist nicht moglich. 
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Die Erfindung geht von der Aufgabe aus, bei einem Lokalisie- 
rungsverf ahren der eingangs genannten Art sowie einer Vor- 
richtung zur Durchfuhrung des Verfahrens die Lokalisierungs- 
genauigkeit beispielsweise mittels Fluorophore f luores zierend 
markierter Tumore zu erhohen sowie eine Auswertung in tiefen 
Gewebeschichten zu ermoglichen und den Rechenaufwand und als 
Folge davon die Rechenzeit drastisch zu senken. 

Die Aufgabe wird erf indungsgemafi fur das Verfahren durch die 
Merkmale des Patentanspruches 1 gelost. Durch das erfindungs- 
gemafie Verfahren kann das Problem der Lokalisierung fluores- 
zierender Objekte in optisch truben Medien schnell gelost 
werden. Weiterhin wird durch die Variation des Anregungsortes 
die Genauigkeit erhoht. 

Die Lichtemission von Gewebeabschnitten, in welchen Fluores- 
zenzmarker angereichert sind, wird durch Einstrahlung von La- 
ser licht geeigneter Wellenlange anger egt . Fluoreszenzlicht 
kann dann an der naheliegenden Hautoberf lache gemessen wer- 
den. 

Um Orte und optische Parameter markierter Gewebeabschnitte zu 
bestimmen wird 

• eine Folge von Fluoreszenzanregungen von der Oberflache, 
beispielsweise von verschiedenen Orten mit verschiedenen 
Modulationsf requenzen (einschliefilich Frequenz Null) , in 
das Gewebe eingestrahlt , und dann 

• das Fluoreszenzlicht mit einer Oder mehreren Anordnungen 
geeigneter, auf der Oberflache verteilter Lichtsensoren 
gemessen, um somit zweidimensionale Messwertverteilungen 
zu erhalten, welche von der Art der Anregung abhangen. 

Es hat sich als vorteilhaf t erwiesen, wenn zur Erzeugung der 
verschiedenen Fluoreszenzeigenschaf ten die Bereiche mit fluo- 
reszierenden Markern (Fluorophore) markiert werden. 
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Die Ortsauf losung wird erhoht, wenn die f luoreszenzanregenden 
Lichtsignalen mit verschiedenen Modulationsf requenzen erzeugt 
und in den Gewebeabschnitt eingestrahlt werden. 

5 In vorteilhaf ter Weise konnen die f luoreszenzanregenden 

Lichtsignalen durch Laserlicht geeigneter Wellenlange einge- 
strahlt werden. 

Zweckmafiigerweise werden die Fiihrungsf elder zunachst normiert 
10 und dann transf ormiert , wobei die Fiihrungsf elder zu orthogo- 
nalen Fiihrungsf eldern transf ormiert werden konnen. Weiterhin 
konnen die orthogonalen Fiihrungsf elder mittels einer Singu- 
'nj*. ' larwertzerlegung aus den Fiihrungsf eldern bestimmt werden. 

15 Erf indungsgemafi konnen die optischen Parameter durch Refe- 
renzmessungen bei nicht-f luoreszenzanregenden Wellenlangen 
mittels Schatzverf ahren bestimmt werden. 

Die Aufgabe wird erf indungsgemaii fur eine Vorrichtung zur 
20 Durchfiihrung des Verfahrens durch die Merkmale des Patentan- 
spruches 8 gelost. 

Ein Messsystem kann beispielsweise 8x8 regular, auf einer 
ebenen Messflache angeordnete Lichtsensoren enthalten. Es 
^5 kann jedoch vorteilhaf t sein, mit mehreren solcher planarer 

Systeme gleichzeitig zu messen. So konnen beispielsweise zwei 
Anordnungen von Lichtsensoren vorgesehen sein, die beidseitig 
des zu untersuchenden Gewebeabschnittes anlegbar sind . Da- 
durch konnen bei Messungen an der weiblichen Brust zwei Mess- 
30 flachen auf entgegengesetzten Seiten der Mamma angelegt wer- 
den. Eine vorteilhaf te Ausgestaltung ist die integrierte An- 
ordnung der Mess flachen in Anpressplatten eines Rontgen- 
Mammographie-Gerates . 

35 Im allgemeinen konnen beliebig gekriimmte bzw. krummbare oder 
flexible Messflachen mit beliebiger Anordnung von Lichtsenso- 
ren zum Einsatz kommen. 
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Die Erfindung ist nachfolgend anhand von in der Zeichnung 
dargestellten Ausf uhrungsbeispielen naher erlautert . Es zei- 
gen: 

5 

Figur 1 Aufnahme zur Erlauterung des Prinzips der optischen 
Fluoreszenz-Bildgebung, 

Figur 2 in einer Ubersichtsdarstellung die wesentlichen 
10 Komponenten eines Gerates zum Lokalisieren und 

Klassif izieren einer fokalen Lasion in einem Gewe- 
^ beabschnitt, 

Figur 3 die wesentlichen Verf ahrensschritte zum Lokalisie- 
15 ren einer fokalen Lasion, 

Figur 4 erf indungsgemafier Applikator mit 8x8 Sensoren sowie 
8 neben der Messflache angeordneten Lichtquellen 
zur Erzeugung des Anregungslichtes , 

20 

Figur 5 ein Doppelsystem mit zwei gegeniiberliegenden Appli- 
katoren, 

Figur 6 eine zweidimensionale Messwertverteilung der Konfi- 
--|p3 guration 1 fur die ersten vier Anregungsorte, 

Figur 7 ein Singularwertspektrum der Konf iguration 1, 

Figur 8 Basismaps der Konf iguration 1, 



30 



Figur 9 Lokalisierungsfunktionen zweier Fluorochrom-mar- 
kierten Lasionen, 



Figur 10 die Lokalisierung unterschiedlich tief er Lasionen 
35 mit einem planaren Messsystem und 
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Figur 11 die Lokalisierung unterschiedlich tiefer Lasionen 
mit zwei gegemiberliegenden planaren Messs ystemen . 

Die Ubersichtsdarstellung in Figur 2 zeigt eine Mess- und 
Auswerteanordnung, womit ein begrenztes, in einem biologi- 
schen Gewebeabschnitt 1 angeordnetes Raumgebiet 2 lokalisiert 
und identif iziert werden kann. Dabei ist vorausgesetzt , dass 
das Raumgebiet 2 eine vom ubrigen Gewebeabschnitt 1 verschie- 
dene Fluoreszenzeigenschaf t besitzt. Diese Voraussetzungen 
sind ausreichend gut erfullt, wenn es sich bei dem biologi- 
schen Gewebeabschnitt 1 urn eine weibliche Brust und bei dem 
begrenzten Raumgebiet 2 urn einen Tumor handelt, dem bei- 
spielsweise ein f luores zierender metabolischer Marker zuge- 
fuhrt wurde, dessen Fluoreszenzeigenschaf ten durch bestimmte 
Enzyme f reigeschaltet werden 

Zur Messanordnung gehort ein Applikator 3 mit einer Vielzahl 
von raumlich verteilt angeordneten Photosensoren sowie dane- 
ben in einer Zeile angeordneten Laserdioden, wie dies noch 
ausfuhrlich beschrieben wird. 

Die Erkennung eines Tumors oder einer anderen markierten Re- 
gion geschieht dann liber Bestrahlung der Region mit Licht der 
Laserdioden in der speziellen Anregungswellenlange des Fluo- 
reszenzf arbstof f es und Detektion des emittierten Lichtes 
durch die Photosensoren in der entsprechenden Emissionswel- 
lenlange des Fluorophors . 

Die Photosensoren und Laserdioden des Applikators 3 sind zum 
einen uber elektrische Verbindungsleitungen 4 mit einer elek- 
trischen Steuervorrichtung 5 und zum anderen tiber elektrische 
Verbindungsleitungen 6 mit einer Messwertauf bereitung 7 ver- 
bunden . 

Mit Hilfe der Steuervorrichtung 5 werden dem biologischen Ge- 
webeabschnitt 1 uber eine Anzahl von K Laserdioden, wobei 
\<K<M ist, Impulse von NIR-Licht zugefuhrt, urn dort eine 
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einen gegebenenf alls vorhanden markierten Tumor zum Fluori- 
zieren anzuregen. 

Zur Lokalisierung und Identif izierung von raumlich begrenzten 
Gebieten 2 wird das von den Gebieten 2 emittierte Licht auf 
der Oberflache des Gewebeabschnittes 1 an M Orten mit Photo- 
sensoren gemessen und einer Auswertung zugeflihrt. 

Die Messwertaufbereitung 7 umfasst z.B. Messverstarker , Fil- 
ter und Analog-Digital-Wandler . Die Messwertaufbereitung 7 
ist mit einem oder mehreren Dateneingangen eines elektroni- 
schen Rechners 8 verbunden. Neben den Messwerten wird dem 
Rechner ein Modell 9 des Gewebeabschnitts 1 zur Verfiigung ge- 
stellt, mit dessen Hilfe die oben erwahnten f luoreszief enden 
Gebiete 2 lokalisiert und identif iziert werden, wie weiter 
unten noch beschrieben ist. Das Ergebnis, z.B. in Form einer 
graphischen Darstellung der Anatomie des Gewebeabschnitts, 
worin der Ort der Lichtquellen und damit der Raumgebiete 2 
markiert ist, erfolgt uber einen Monitor 10. Da die Berech- 
nung unter anderem von dem Modell 9 und dem Ort der Beleuch- 
tung bestimmt ist, ist eine iibergeordnete Eingabe und Steue- 
rung 11 vorgesehen, mit der die Anzahl und der Ort der Photo- 
sensoren sowie der Laserdioden, der Wert der Frequenz und das 
Modell vorgegeben werden. 

Das Lokalisierungsverf ahren wird beispielhaft anhand von Fi- 
gur 3 erlautert. Zunachst werden seine Eingabegrofien, d.h. 
die Mess- und die Modelldaten, und dann die Rechenschritte 
des Verfahrens erlautert. 

Die Eingabegroflen fur das Lokalisierungsverf ahren sind pro 
Messf lache 

a) Eine MxN Datenmatrix D mit Messwerten (Bezugszeichen 21), 
welche von den M Sensororten r Sm , (/w = l ) ... 5 M)und den N Anre- 

gungsparametern ( N x Anregungsorte (n x = \... 7 N l ) und/oder 

N 2 Anregungsmodulationsf requenzen / (n 2 = 1,...,N 2 ) , wobei 
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N = N l +N 2 ), abhangt, und welche sich ggf s . durch Nachver- 
arbeitung aus den eigentlichen Messdaten ergeben kann, 

Die M-dimensionalen Spaltenvektoren der Datenmatrix las sen 
sich entsprechend der Anordnung der Sensoren auf der Mess- 
flache ref ormatieren . Eine graphische Darstellung des re- 
formierten Spaltenvektors visualisiert die Messwertvertei- 
lung liber der betrachteten Messflache zu gegebener Anre- 
gungsart. Im Falle oben erwahnter 8x8 Sensorverteilung 
wird der 64-dimensionale Spaltenvektor in eine 8x8-Matrix 
umgef ormt . 

- von optischen Parametern wie den Absorptions- und Streuko- 
effizienten ju a ,ju s des die Lasion(en) umgebenden Mediums. 

b) Ein Satz von K Fuhrungsf eldern oder Leadfields 
L k (f m9 n m9 r i9 ^ a9 /i M ) 9 (k = l 9 ... 3 K) i beispielsweise Multipol- 

Leadfields, welche in Figur 3 mit dem Bezugszeichen 22 ge- 
kennzeichnet sind und welche ihrerseits abhangen 

- vom Modell des optischen Mediums des Untersuchungsgebiets 
1/ 

- vom Messsystem, beispielsweise Ort f m und/oder Normalen- 
vektor n m des m-ten Sensors, 

- vom Ort f f des T-ten anregbaren Fluorochroms, 

- von der Art der Messung ( Frequenzmodulation ja/nein) und 

- von optischen Parametern wie den Absorptions- und Streuko- 
effizienten ju a >M s des die Lasion(en) umgebenden Mediums 

Die auf einen Sensor auf tref f enden Photonen werden in elekt- 
rische Signale umgewandelt und dann der weiteren Auswertung 
zugefiihrt. Im Falle f requenzmodulierter Anregung werden 
Lichtintensitat und Phasenverschiebungen gegenuber der Ein- 
gangswelle gemessen. Beide reellen Messwerte konnen zu einem 
komplexen Messwert zusammengef asst werden. Die Datenmatrix 
ist dann - im mathematischen Sinne - komplex. Im Folgenden 
wird vom allgemeinen Fall einer komplexen Datenmatrix ausge- 
gangen. 
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Es kann erf orderlich sein, nachverarbeitete Messdaten dem Lo- 
kalisierungsalgorithmus zuzufiihren. Beispielsweise werden 
durch Abschneiden von Randdaten Randartef akte eliminiert. Sie 
konnten eine nicht-existente Abhangigkeit von der Modulati- 
5 onsf requenz- bzw. dem Anregungsort vortauschen . 

Die Datenmatrix kann sich auch aus einer Linearkombination 
von mindestens zwei Datensatzen ergeben. Beispielsweise kann 
die Differenz eines Datensatzes mit Fluoreszenzsignalen und 
10 eines raumlich benachbarten Datensatzes ohne Fluoreszenzsig- 

nal betrachtet werden. Es ist zu erwarten, dass mogliche Bei- 
'-0 trage von Hintergrundanregungen in den Dif f erenzdaten deut- 
lich reduziert, wenn nicht ganzlich eliminiert sind. 

15 Fiihrungsf elder oder sogenannte Leadfields sind aus dem Bio- 
elektroraagnetismus bekannte Grofien . Sie geben die Messwert- 
verteilung einer Einheitssignalquelle an, die mit einem gege- 
benen Messsystem aufgenommen werden kann, 

20 Leadfields, welche die mit einem Messsystem oder mehreren 
Messsystemen auf nehmbare Lichtintensit&t auf Grund optisch 
angeregter, fokaler und mit Fluorochromen markierten Lasionen 
beschreiben, sind als Eingangsgrofien fur das Verfahren zur 
Lokalisierung solcher fokaler Lasionen geeignet. 

Beispielhaft wird im Ausf uhrungsbeispiel lediglich ein Lead- 
field verwendet . Es beschreibt die mit einem gegebenen Mess- 
system rnessbare Lichtintensitat einer punktf ormigen Licht- 
quelle. Raumlich ausgedehnte Fluoreszenzquellen konnen ent- 

30 sprechend den in B. Scholz, "Towards Virtual Electrical 

Breast Biopsy: Space-Frequency MUSIC for Trans-Admittance Da- 
ta", IEEE Trans. Med. Imag., Vol. 21, No. 6, pp. 588-595, 
2002, behandelten ausgedehnten elektrisch-polarisierten Lasi- 
onsgebiete ebenfalls durch Multipol-Leadf ields erfasst wer- 

35 den. Im folgenden wird angenommen, dass beispielhaft Multi- 
pol-Leadf ields , d . h . ein Satz von mehreren Leadfields , zur 
Verfiigung stehen. 
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Fur die weiteren Schritte ist es hilfreich, die Werte des k- 
ten Leadfields L k (k = \ y ... 7 K) an den M Messorten zu einem M-di- 

mensionalen Vektor im Datenraum ( symbolisiert durch den Un- 
5 terstrich unter L) zusammenzuf assen. 

0) L k (r)^(L k (r,r l ) 9 ... 9 L k (r 9 r M )J mit k = l 9 ...,K 



Hierbei ist f der Schwerpunktsort der Lasion. Der Obersicht- 
10 lichkeit halber ist in Gleichung (1) die Abhangigkeit von den 

optischen Parametern des die Lasion (en) umgebenden Mediums 
^ nicht aufgefuhrt. 

Die optischen Parameter, welche wie oben erwahnt in die Lead- 
15 fields eingehen, konnen durch Ref erenzmessungen bei nicht- 
f luores zenzanregenden Wellenlangen mitt els Schatzverf ahren 
bestimmt werden. 



Die Signalverarbeitung des Verfahrens besteht pro Messflache 
2 0 aus 

1 . der Singularwertzerlegung der Datenmatrix D (Bezugszeichen 
23 in Figur 3) , 

2. der Analyse der Singularwertzerlegung (Bezugszeichen 24 in 
^ Figur 3) und 

,jw5 3. dem eigentlichen Lokalisierungsverf ahren (Bezugszeichen 25 
in Figur 3) . 



Die Singularwertzerlegung 28 einer Matrix ist eine aus G. Go- 
lub, Ch. Van Loan, Matrix Computations, 3 rd edition, J. Hop- 
30 kins University Press, 1996, Seite 70 ff., bekannte mathema- 
tische Methode. Sie lautet fur obige Datenmatrix 

(2) D = USV H . 



35 Hierbei bezeichnet 

U eine nur von den Indizes der Sensororte abhangi- 
ge, unitare MxM Matrix, 
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S die MxN Singularwert-Matrix mit min(M,N) reellen 
Singularwerten in der Diagonalen und sonst ver- 
schwindenden Elementen und 

V eine nur von den Anregungsorts- bzw. Frequenzin- 
dizes abhangige, unitare NxN Matrix und 

H die hermitesche Konjugation der betreffenden Mat- 
rix . 

Die Singularwerte sind entsprechend ihrer abnehmenden numeri- 
schen Grofie geordnet, d.h. es gilt 

C*) s } > s 2 > ...> s^m^ 

Bezeichnet man mit H q ,^ q die g-ten Spaltenvektoren der Matri- 

zes U und V, dann zeigt die alternative tensorielle Schreib- 
weise (® bezeichnet das Tensorprodukt ) 

(4) D = ^^s q u q ^f q 

deutlich, dass der gr-te Singularwert ausschlieftlich mit den 
qr-ten Spaltenvektoren von U und V verknupft ist. Der einfache 
und der doppelte Unterstrich bei u und v sollen andeuten, 
dass es sich urn einen M- bzw. N-dimensionalen Vektor handelt . 

Die M Indizes der Spaltenvektoren u q entsprechen den fort- 

laufend nummerierten Indizes der Messsensoren . Demzufolge 
konnen diese Spaltenvektoren - wie oben erwahnt - entspre- 
chend der Anordnung der Messsensoren in Matrizes umgeformt 
und als zweidimensionale Messwertverteilungen dargestellt 
werden. Diese Spaltenvektoren sind anregungs- bzw. frequenz- 
unabhangige orthonormierte Basisvektoren im M-dimensionalen 
Datenraum und werden hier als Basis- bzw. Eigenmaps bezeich- 
net . 

Zur Singularwertanalyse wird die Zahl Q dom der signif ikanten 
Singularwerte ermittelt, die die Zahl der sich bzgl. der An- 
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regungsart linear unabhangig verhaltenden Fluoreszenzquellen 
angibt . 

Eine punktf ormige Inhomogenitat im ansonsten homogenen opti- 
5 schen Medium erzeugt beispielsweise ein Singularwertspektrum 
mit einem signif ikanten Singularwert (Q dom = l) - 

Die zugehorigen Spaltenvektoren u q werden als Basisvektoren 

eines - f requenzunabhangigen - Odo^-dimensionalen Signalraumes 
10 im M-dimensionalen Datenraum betrachtet. Die restlichen M- 

Qdom Spaltenvektoren sind dann die Basisvektoren des orthogo- 
Vg^' nalen Signalraumes. 

Das Aufsuchen von Fuorochrom-markierten Lasionen, die Lokali- 
15 sierung, entspricht der Suche von Orten bzw. Schwerpunktsor- 
ten angeregter Signalquellen . Diese Suche mittels eines Com- 
puters verlangt die Diskretisierung des angenommenen Modell- 
mediums, welche die zu untersuchende Korperregion mathema- 
tisch nachbilden soli. 

20 

Eine Suchstrategie besteht darin, mit den anregungs- und fre- 
quenzunabhangigen Leadfields an jedem Rasterort anregungs- 
und f requenzunabhangige Modelldaten und/oder einen Modellda- 
tenraum zu erzeugen und diese und/oder diesen mit dem aus den 
;-^p5 Messdaten gewonnenen anregungs- und f requenzunabhangigen Sig- 
nalraum zu vergleichen. Vergleichsmafie lassen sich so defi- 
nieren, dass sie den Grad der „Ubereinstimmung* zwischen Sig- 
nalraum und Modelldaten / Modelldatenraum anzeigen. Orte, an 
denen das Maft ein lokales Maximum annimmt, werden als Orte 
30 tatsachlicher Signalquellen angesehen. 

Eine alternative, zweite Suchstrategie besteht im Vergleich 
zwischen dem orthogonalen Signalraum - in der alteren Litera- 
tur auch Rauschraum genannt - und den Modelldaten bzw. dem 
35 Modelldatenraum. Vergleichsmafie lassen sich dann so definie- 
ren, dass sie den Grad der „Nicht-Ubereinstimmung* zwischen 
dem orthogonalen Signalraum und Modelldaten / Modelldatenraum 
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anzeigen. Orte, an denen das MaB ein lokales Minimum annimmt, 
werden als Orte tatsachlicher Signalquellen angesehen . 

Die Modelldaten sind durch die Leadfields gegeben: entweder 
5 werden sie direkt oder nachverarbeitet genutzt. 

Ein einzelnes Leadfield stellt ein Modelldatensatz dar, der 
eine spezielle Eigenschaft der Signalquelle widerspiegelt . 
Beispielsweise beschreibt das Leadfield einer punktf ormigen 
10 Fluoreszenzquelle, die messbare Lichtintensitat bei isotroper 
Lichtemission durch diese Quelle, 

Die Gesamtheit der betrachteten Leadfields (Anzahl: K) defi- 
niert auf Grund ihrer linearen Unabhangigkeit einen K-dimen- 

15 sionalen Modelldatenraum. Die Leadfields sind i.a. nicht- 

orthogonale Basisvektoren dieses Modelldatenraumes . Orthogo- 
nale Basisvektoren lassen sich durch geeignete Orthogonali- 
sierungsverfahren, d.h. durch Nachverarbeitung der Lead- 
fields, gewinnen. Sie verandern den Model ldatenraum nicht. 

20 Jedoch ergeben sich mit den neuen Basisvektoren neue einzelne 
Modelldatensatze (s. oben) . Diese Basisvektoren konnen zu- 
satzlich normiert werden. Dies stellt sicher, dass Leadfields 
mit unterschiedlichem Abstandsverhalten in gleicher Weise zur 
Lokalisierung beitragen konnen. Zudem hat es den Vorteil, 

0E5 physikalisch dimensionslose Grolien zu betrachten. 

Eine vorteilhafte Leadf ield-Nachverarbeitung besteht bei- 
spielsweise darin, die K Leadfields L k [k=l,...,K) aus Glei- 

chung (1) zu normieren ( Verarbeitungsschritt 27) . Dabei wer- 
30 den jeweils die einzelnen Fuhrungsf elder auf ihre Norm bezo- 
gen, so dass sich die normierten Fuhrungsf elder L^ n ^ wie 
folgt ergeben: 



(4a) 




T (n) _ L k 
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Mittels beispielsweise einer Singularwertzerlegung der MxK 
Leadf ields-Matrix L orthogonalisierte Leadfields zu gewinnen. 
Die Normierung ist durch den Index (n) angezeigt. 

(5) L w = ^,...,Lf)=U L S L Vl 

Der Ubersichtlichkeit halber wurden die Argumente der Lead- 
fields, die Ortsvektoren des Quellortes, weggelassen. Die 
ersten K Spaltenvektoren U(r) Lk ,(k = l,...,K) der Matrix U L sind 

die gesuchten quel lortsabhangi gen orthonormierten Leadfields. 
Im Falle eines einzigen Leadfields entfallt die Singularwert- 
zerlegung aus Gleichung (5) . 

Beispiele fur Vergleichsmaiie zwischen einem Modelldatensatz 
bzw. dem Modelldatenraum und dem Signal- bzw. dem orthogona- 
len Signalraum sind aus anderen biomedizinischen Anwendungen, 
Analyse biomagnetischer Daten oder Analyse elektrischer 
Trans-Admittanzdaten, bekannt. Solche Verfahren sind Projek- 
tionsverf ahren und Winkelabstandsverfahren. 

Mit Hilfe von Pro j ektionsmatrizes werden einzelne Modellda- 
tensatze bzw. der Modelldatenraum entweder auf den Signalraum 
bzw. auf den orthogonalen Signalraum projiziert und fiir jeden 
Rasterort der entsprechende Proj ektionswert bestimmt. 

Basierend auf dem in G. Golub et al. Seite 584 ff. angegebe- 
nen Algorithmus zur Berechnung von Winkeln zwischen zwei Un- 
terraumen, dem sogenannten Winkelverf ahren, werden Suchort 
fur Suchort die Winkel zwischen dem Signalraum bzw. dem or- 
thogonalen Signalraum und einzelnen Modelldatensat zen bzw. 
dem Modelldatemraum berechnet. Hier gibt ein kleiner Winkel, 
also ein kleiner Wert des VergleichsmaBes, zwischen bei- 
spielsweise dem Signalraum und dem Modelldatenraum eine grofie 
„Obereinstimmung XN an. Eine Transformation des Vergleichsmalies 
in der Form 90 Grad - berechneter Winkel ergibt dann wieder 
Maxima der Vergleichsfunktion an den Orten tatsachlicher Sig- 
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nalquellen. Mutatis mutandis lasst sich das Gesagte auf Win- 
kelvergleichsmaJSe zwischen anderen Unterraumen iibertragen. 

An jedem Ort f des diskretisierten optischen Model lmediums 
5 wird uberpruft, wie grofi der Abstand zwischen den orthogona- 
lisierten Leadfields U(r) Lk und dem Signalraum ist. Ein ge- 

eignetes Mali ist die Funktion 
(6) F k (r) = (ZlT c ^<-^. k ] 2 

10 

Die Ausgangsgleichung von (6) ist die im Quadratmittelsinne 
/ )|| zu betrachtende Gleichung 

(?) ZS~c,t/,=C/ i>k k = l,...,K. 

15 

Wird die Losung fur die Koef f izienten c t in das Bewertungsmafi 
eingesetzt, dann folgt 

(8) ^O0=i-z£r<kf m\*? • 

20 

Dieses Mafi entspricht einer Projektion des betrachteten Lead- 
fields auf den orthogonalen Signalraum. Unter Benutzung der 
auf den orthogonalen Signalraum pro j izierenden Projektions- 
matrix 

25 

(9) P a ,=l-'Zu l ®u? 
ergibt sich 

30 (10) F k (r) = \P os U(r) LJ \* 

Die tatsachliche Lokalisierungsfunktion F ist der Minimalwert 
der Abstande F k . Sie ist definiert durch 

35 (11) Ffr) = ttm{F k (r)) 
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Die lokalen Minima der Lokalisierungsfunktion werden entspre- 
chend ihrer Zahlenwerte monoton aufsteigend geordnet. Die Or- 
te, welche den ersten Q dom lokalen Minima zuzuordnen sind, 

werden als Orte von Signalgeneratoren angesehen. 

Im Falle mehrerer (M^ s ) Messflachen werden die obengenannten 

Rechenschritte fur die Daten jeder Messflache separat ausge- 
ftihrt. Pro Messflache ergibt sich dann eine Zielfunktion ge- 
mafi Gleichung (11) . Aus diesen Einzelzielfunktionen lasst 
sich eine Gesamtzielfunktion F igesamt) gemafi 



definieren. Hierbei ist F iM) die Zielfunktion der //-ten Mess- 
flache . 

Die lokalen Minima der Gesamt-Lokalisierungsfunktion werden 
wie oben entsprechend ihrer Zahlenwerte monoton aufsteigend 
geordnet. Die Orte, welche den ersten Q dom lokalen Minima zu- 
zuordnen sind, werden als Orte von Signalquellen angesehen. 
Das Ausfuhrungsbeispiel bestatigt die Erwartung, dass bei 
mehreren nichttrivial angeordneten Messflachen die lokalen 
Minima der Einzelzielfunktionen deutlicher ausgepragt sind 
und damit das Lokalisierungsergebnis sicherer machen. 

Die Auf fiihrungsbeispiele wurden mit planaren, in dem Applika- 
tor 3 angeordneten Messsystemen gewonnen, welche 8x8 regular 
angeordnete Photosensoren 31 enthielten, wie sie in Figur 4 
schematisch dargestellt sind. Die Sensoren 31 wurden punkt- 
formig angenommen. Ihr Abstand langs einer Richtung betrug 8 
mm, so dass sich eine Messf eldf lache 56x56 mm 2 ergibt . Die 
Orte, an denen sich 8 Laserdioden 32 befinden, welche das 
Fluoreszenz anregende NIR-Licht in die Korperregion einstrah- 
len, konnen beispielsweise neben der Messflache angeordnet 
sein. Die Anregung kann, muss aber nicht, f requenzmoduliert 
sein. Eine solche Messanordnung kann mit der Hand uber einen 



(12) 
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interessierenden Gewebeabschnitt 1 gefuhrt werden. Die Laser- 
dioden 32 senden Anregungsstrahlen 33 aus, die auf das fluo- 
reszierende Raumgebiet 2 treffen. Die Fluoreszensstrahlen 34 
werden von den Photosensoren 31 erf asst. 

In Figur 5 ist ein Doppelssytem eines Applikators 3 mit zwei 
gegenuberliegenden planare Messflachen gleicher Dimensionie- 
rung (8x8 Sensoren) dargestellt. Sie konnen beispielsweise in 
den Anpressplatten eines Rontgenmammographiegerates integ- 
riert sein. Die Fluoreszenzanregung erfolgt an 8 Anregungsor- 
ten, welche sich neben der Messflache (z=0) des oberen Appli- 
kators 3 befinden. Der Abstand der beiden Applikatoren 3 be- 
tragt 64 mm. 

Als optische Gewebemodelle werden bei den vorliegenden Bei- 
spielen folgende Modelle herangezogen : 

A) Das einfachste Modell ist ein berandungsloses Gebiet mit 
punktf ormigen f luoreszierenden Objekten, welches ansonsten 
optisch homogen (konstante optische Parameter wie Absorp- 
tions- und Streukoef f izient ) ist. 

B) Als zweites Modell wurde ein optisch inhomogener Quader 
Gebiet mit punktf ormigen f luoreszierenden Objekten be- 
trachtet. Es wurde angenommen, dass Absorptions- und 
Streukoef f izient lokal um 100% variieren konnen. Die Figur 
9 zeigt die Lokalisierungsf unktionen von 32 mm und 48 mm 
tie fen, Fluorochrom-markierten Lasionen . Der Absorpti- 
onskonstrastunterschied des umgebenden Gewebes betragt 
100% (Bild in Bild) . 

Der Simulation der Daten werden folgende Konf igurationen 
zugrunde gelegt: 
Konfi gura tion 1 

Mess-/Anregungssystem: siehe Figur 4, die Anregung ist 

nicht f requenzmoduliert 
Gewebemodell : inhomogener Quader (5.2.B) 
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Fluoreszenzquelle : 



Daten: 
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Ort bei (x, y, z ) = (28 , 28 , 32 ) mm, 
d.h. zentrale Lage unterhalb der 
Messflache in einer Tiefe von 32 
mm (Koordinatensystem s . Figur 4) 
siehe Figur 6 



10 



15 



20 



Kon figuration 2 
Mess-/Anregungs system: 

Gewebemodell : 
Fluoreszenzquelle : 



Konfiguration 3a A 4a und 5a 
Mess-/Anregungs system: 



Gewebemodell 



Einzel-Fluoreszenzquellen : 



siehe Figur 4, die Anregung ist 
nicht f requenzmoduliert 
inhomogener Quader (5 . 2 .B) 
Ort bei (x, y, z) = (28, 28, 48 ) mm, 
d.h. zentrale Lage unterhalb der 
Messflache in einer Tiefe von 48 
mm (Koordinatensystem s. Figur 4) 



Einzelmesssysterrw siehe Figur 4, 
die Anregung ist nicht frequenz- 
moduliert 

homogenes, unberandetes Medium 
(5.2 .A) 

Orte bei (x, y, z) = (28, 28, 16) mm, 
(28,28,32) mm, (28,28,48) mm d.h, 
zentrale Lagen unterhalb der 
Messflache in Tiefen von 16 mm, 
32 mm und 48 mm (Koordinatensys- 
tem s . Figur 4 ) 



Konfiguration 3b, 4b und 5b 
30 Mess-/Anregungs system: 



Gewebemodell : 



35 Einzel-Fluoreszenzquellen ; 



Doppelmessystem, siehe Figur 5, 
die Anregung ist nicht frequenz- 
moduliert 

homogenes, unberandetes Medium 
(5.2 .A) 

Orte bei (x, y, z) = (28, 28, 16) mm, 
(28,28,32) mm, (28,28,48) mm d.h, 
zentrale Lagen unterhalb der 



Copied from 10691407 on 06/17/2005 

::::::::::: :::::::::::: . 



200219266 



20 

Messflache in Tiefen von 16 mm, 
32 mm und 48 mm (Koordinatensys- 
tem s. Figur 5) 

Das in Figur 7 dargestellte Singularwertspektrum der Daten 
der Konf iguration 1 aufgrund der Singularwertzerlegung 23 
weist erwartungsgemafi - entsprechend der Zahl der vorhandenen 
Fluoreszenzquellen - einen numerisch dominanten Singularwert 
auf. Die restlichen Singularwerte geben Rauschen, in diesem 
Fall numerisches Rauschen, wieder. 

Die zugehorigen Basismaps bzw. Eigenmaps zeigt die Figur 8. 
Entsprechend dem einzigen numerisch dominanten Singularwert 
gibt es eine strukturierte Basismap. Sie definiert den hier 
eindimensionalen Signalraum des (hier 64-dimensionalen) Da- 
tenraumes . 

Die oben definierten Zielfunktionen zur Lokalisierung 25, 
d.h. die Lokalisierungsfunktionen, der Konf igurationen 1 und 
2 sind in Figur 9 dargestellt. 

Der Einfluss einer zweiten gegenuberliegenden Messflache auf 
die Lokalisierung wird an Hand der Zielfunktionen der Konfi- 
gurationen ly a/b, 4 a/b und 5 a/b gezeigt, wobei Figur 10 die 
Lokalisierung unterschiedlich tiefer Lasionen mit einem pla- 
naren Messsystem und Figur 11 die Lokalisierung unterschied- 
lich tiefer Lasionen mit zwei gegenuberliegenden planaren 
Messsystemen zeigen. Die Lagen der Messsonden sind durch di- 
cke Striche an dem linken Rand bzw. an beiden seitlichen Ran- 
dern markiert. Gegenuber Figur 10 ist eine deutlichere Aus- 
pragung der Minima sichtbar, wobei zu beachten ist, dass die 
Skala gemafi Figur 11 gegenuber der von Figur 10 unterschied- 
lich ist. 

Durch das erf indungsgemaiie Verf ahren kann das Problem der Lo- 
kalisierung f luoreszierender Objekte in optisch truben Medien 
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schnell gelost werden. Weiterhin wird durch die Variation de 
Anregungsortes die Genauigkeit erh5ht. 

Dieses Lokalisationsverf ahren zeichnet sich dadurch aus, das 
es 

• in Echtzeit arbeitet, 

• patient enunabhangig ist und 

• robust gegen Schatzung optischer Parameter ist. 
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Patentanspruche 

1. Verfahren zur Lokalisierung von Bereichen (2) in einem 
biologischen Gewebeabschnitt (1), die zumindest wahrend der 
Untersuchung eine vom Gewebeabschnitt (1) verschiedene Fluo- 
reszenzeigenschaf t aufweisen, aufgrund derer bei einer Be- 
strahlung mit Licht (33) einer ersten Wellenlange Licht (34) 
einer anderen Wellenlange emittiert wird, mit den Schritten: 

a) Anlegen einer Folge von f luores zenzanregenden Lichtsigna- 
len an unterschiedlichen Orten auf dem Gewebeabschnitt 
(1), 

b) Messen von Fluoreszenzlicht (21) an mehreren Messorten auf 
einer Oberflache des Gewebeabschnittes (1), die sich auf- 
grund der Lichtsignale dort einstellen, 

c) Bestimmen von f requenzunabhangigen Signalanteilen in den 
Antwortsignalen und Weiterverarbeitung (23, 24) der fre- 
quenzunabhangigen Signalanteile zu Eingabewerten eines Lo- 
kalisationsschritts (25) , 

d) Modellieren des Gewebeabschnitts (1) und Bestimmen eines 
Satzes von Fuhrungsf eldern (22) , 

e) Transf ormierung (28) der Fuhrungsf elder (22), so dass im 
Lokalisierungsschritt (25) die f requenzunabhangigen Sig- 
nalanteile mit den transf ormierten Fuhrungsf eldern vergli- 
chen werden und dass der Ort (2 6) der transf ormierten Fuh- 
rungsf elder, die die f requenzunabhangigen Signalanteile am 
besten wiedergeben, als Ort des zu lokalisierenden Be- 
reichs (2) ausgegeben wird. 

2. Verfahren nach Anspruch 1, dadurch g e - 
kennzeichnet, dass zur Erzeugung der verschie- 
denen Fluoreszenzeigenschaf ten die Bereiche (2) mit fluores- 
zierenden Markern (Fluorophore) markiert werden. 

3. Verfahren nach Anspruch 1 oder 2, dadurch 
gekennzeichnet , dass die f luoreszenzanregen- 
den Lichtsignalen mit verschiedenen Modulationsf requenzen er- 
zeugt und in den Gewebeabschnitt (1) eingestrahlt werden. 
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4. Verfahren nach einem der Ansprtiche 1 bis 3, d a - 
durch gekennzeichnet, dass die fluo- 
reszenzanregenden Lichtsignalen durch Laserlicht (33) geeig- 
neter Wellenlange eingestrahlt werden. 

5. Verfahren nach einem der Ansprtiche 1 bis 4, d a - 
durch gekennzeichnet, dass die Fiih- 
rungsf elder (22) zunachst normiert (27) und dann transfor- 
miert (28) werden. 

6. Verfahren nach einem der Ansprtiche 1 bis 5, d a - 
durch gekennzeichnet, dass die Ftih- 
rungsfelder zu orthogonalen Ftihrungsf eldern transf ormiert 
(28) werden. 

7. Verfahren nach Anspruch 6, dadurch ge- 
kennzeichnet , dass die orthogonalen Ftihrungs- 
f elder mittels einer Singularwertzerlegung (28) aus den Ftih- 
rungsf eldern (22) bestimmt werden. 

8. Verfahren nach Anspruch 7, dadurch ge- 
kennzeichnet , dass die optischen Parameter 
durch Ref erenzmessungen bei nicht-f luoreszenzanregenden Wei- 
lenlangen mittels Schat zverf ahren bestimmt werden. 

9. Vorrichtung zur Durchflihrung des Verfahrens nach einem 
der Ansprtiche 1 bis 8, gekennzeichnet 
durch, dass wenigstens eine Anordnung (3) von auf 
der Oberflache des Gewebeabschnitts (1) verteilter Lichtsen- 
soren (31) zur Messung des von dem f luoreszierend markierten 
Bereich (2) emittierten Fluoreszenzlichtes (34) vorgesehen 
ist, neben denen Laserdioden (32) zur Erzeugung von den fluo- 
reszierend markierten Bereich (2) anregendem Licht (33), um 
somit von der Art der Anregung abhangende zweidimensionale 
Messwertverteilungen zu erhalten. 
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10. Vorrichtung nach Anspruch 9, dadurch ge- 
kennzeichnet, dass zwei Anordnungen (3) von 
Lichtsensoren (31) vorgesehen sind, die beidseitig des zu un- 
tersuchenden Gewebeabschnittes (1) anlegbar sind. 

11. Vorrichtung nach Anspruch 8 oder 10, dadurch 
gekennzeichnet, dass die Anordnungen (3) in 
Anpres splat ten eines Rontgen-Mammographie-Gerates integriert 
sind. 

12. Vorrichtung nach einem der Anspruche 8 bis 11, da- 
durch gekennzeichnet, dass die An- 
ordnungen (3) flexibel ausgefuhrt sind. 

13. Vorrichtung nach einem der Anspruche 8 bis 12, da- 
durch gekennzeichnet, dass die An- 
ordnungen (3) gekrummt ausgefuhrt sind. 
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Zusammenf as sung 

Verfahren und Vorrichtung zur Lokalisierung von Licht emit- 
tierenden Bereichen 

Die Erfindung betrifft ein Verfahren sowie eine Vorrichtung 
zur Durchfuhrung des Verfahrens zur Lokalisierung von Berei- 
chen (2) in einem biologischen Gewebeabschnitt (1), die zu- 
mindest wahrend der Untersuchung eine voia Gewebeabschnitt (1) 
verschiedene Fluoreszenzeigenschaf t aufweisen, aufgrund derer 
bei einer Bestrahlung mit Licht (33) einer ersten Wellenlange 
Licht (34) einer anderen Wellenlange emittiert wird, mit den 
Schritten: 

a) Anlegen einer Folge von f luoreszenzanregenden Lichtsigna- 
len (33) an unterschiedlichen Orten auf dem Gewebeab- 
schnitt ( 1 ) , 

b) Messen (21) von Fluoreszenzlicht (34) an mehreren Messor- 
ten auf einer Oberflache des Gewebeabschnittes (1), die 
sich aufgrund der Lichtsignale dort einstellen, 

c) Bestimmen von f requenzunabhangigen Signalanteilen in den 
Antwortsignalen und Weiterverarbeitung (23, 24) der fre- 
quenzunabhangigen Signalanteile zu Eingabewerten eines Lo- 
kalisationsschritts (25) , 

d) Modellieren des Gewebeabschnitts (1) und Bestimmen eines 
Satzes von Ftihrungsf eldern (22), 

e) Transf ormierung (28) der Ftihrungsf elder (22), so dass im 
Lokalisierungsschritt (25) die f requenzunabhangigen Sig- 
nalanteile mit den transf ormierten Ftihrungsf eldern vergli- 
chen werden und dass der Ort (2 6) der transf ormierten Fuh- 
rungsfelder, die die f requenzunabhangigen Signalanteile am 
besten wiedergeben, als Ort des zu lokalisierenden Be- 
reichs (2) ausgegeben wird. 

Figur 2 
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FIG 2 
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FIG 6 



FIG 7 
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• FIG 9 



Hand-held Probe 
Localization of Fluorescence Spots, Depths: 32 mm, 48 mm 
Inhomgeneous Cuboid Medium 
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FIG 10 




z / mm (^Distance from Upper Probe at z=0) 
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